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Observation of the growth of a magnetic vortex in the transition layer of a mildly 

relativistic oblique plasma shock 
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A 2D particle simulation models the collision of two electron-ion plasma clouds along a quasi- 
parallel magnetic field. The collision speed is 0.9c and the density ratio 10. A current sheet forms at 
the front of the dense cloud, in which the electrons and the magnetic field reach energy equipartition 
with the ions. A structure composed of a solenoidal and a toroidal magnetic field grows in this sheet. 
It resembles that in the cross-section of the torus of a force- free spheromak, which may provide the 
coherent magnetic fields in gamma-ray burst (GRB) jets needed for their prompt emissions. 
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The ultrarelativistic jet of a Gamma Ray burst (GRB) 
gives rise to powerful eruptions of electromagnetic ra- 
diation, which can be detected across cosmological dis- 
tances. The non- stationary plasma acceleration at the 
source of the GRB jet results in a nonuniform plasma 
density and flow speed within the jet. Plasma clouds will 
thus collide at a mildly relativistic speed, triggering the 
formation of radiative shocks at the boundaries between 
individual clouds. The ensemble of these shocks is, ac- 
cording to the internal shock model [l|, l2| , the source of 
the "prompt" GRB emissions. 

Plasma processes and instabilities in the shock transi- 
tion layer provide the relativistic electrons and the mag- 
netic fields, which result in the emission of the electro- 
magnetic radiation The filamentation instability 
and the formation of relativistic shocks has been widely 
investigated numerically with the help of particle-in-cell 
(PIC) simulations for electron-positron plasmas [6|-l8| and 
for electron- ion plasmas |9l-[llj, which were all initially 
unmagnetized. However, the magnetic fields driven by 
the filamentation instability in an initially unmagnetized 
plasma tend to have too short a lifetime , and too small 
a coherence length to explain the prompt GRB emissions 

The linear polarization of the GRB emissions sug- 
gests the presence of large-scale magnetic fields within 
the jets [13]. This serves as a motivation to examine 
further with PIC simulations the shock formation and 
evolution in the presence of an ambient magnetic field. 
Flow- aligned [14] and oblique [l5l-[l7| guiding magnetic 
fields have been introduced in previous PIC simulations 
of electron-positron and electron-ion shocks. In particu- 
lar the oblique electron-ion shocks are formidable ampli- 
fiers of the magnetic fields and accelerators of the elec- 
trons. However, the studies in the Refs. 

Ella, 

that 

could be representative for a magnetized internal GRB 
shock, employed a one-dimensional simulation geometry, 
which suppresses the beam filamentation and, as we show 
here, the formation of current vortices. 

A PIC simulation code solves the Maxwell's equa- 
tions and the relativistic Lorentz force equation for the 




FIG. 1: Initial conditions: The simulation box is subdivided 
into two halves, one for each cloud. The clouds collide with 
the speed v c = 0.9c along x. The magnetic field Bo is oriented 
in the x—z plane. The convection electric field Eo points along 
y and changes its sign at x = 0. 



computational particles [18) . We perform here a two- 
dimensional simulation study with the widely used PSC 
code [19]. Figure Q] shows the simulation setup. Two 
plasma clouds, each consisting of ions and electrons with 
the mass ratio rrii/m e = 250, collide at the position 
x = 0. The initial electron and ion number densities 
of the dense cloud are both equal to m and those of the 
tenuous cloud are both equal to ri2 = ni/10. The velocity 
vectors of both clouds are antiparallel and aligned with x. 
The modulus of each cloud in the box frame gives the 
collision speed v c = 2^/(1 + v 2 /c 2 ) = 0.9c. The dense 
cloud propagates to increasing values of x. The temper- 
ature of all 4 plasma species is 131 keV, which yields an 

1/2 

electron thermal speed v t = (ksT '/m e ) 1 , which is close 
to Vb. The cloud density ratio and the collision speed are 
both compatible with the internal shock model [2j. 

The plasma frequency u P 2 — (nie 2 /ra^eo) 1 ^ 2 of the 
ions of the dense cloud normalizes all quantities, which 
have the subscript p for MKS units. The time t = uj P 2 t p , 
the position (x,y) = (x p ,y p )/\i with A$ = c/u P 2, the 
velocity v = v p /c, the electric and magnetic fields 
E = e~Ep/(rriiCUJp2) and B = e3 p /(rriiUJ p 2). The mo- 
menta p = (p x ,Py,Pz) are normalized as p = p^/m^c 
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(ions) and p = p p /m e c (electrons). The initial B 0p gives 
a uj ce = e|Bo p |/m e that equals the electron plasma fre- 
quency uj p i = 250 1 / 2 cj P 2 of the dense cloud. We obtain 
| B 1 = 250 1 / 2 with B 0x = 10B 0z . The modulus of the 
convection electric field is \Eo y \ = v^B^z/c. The sim- 
ulation box size L x x L y = 656A^ x 6A^ is resolved by 
2.8 • 10 4 x 256 rectangular grid cells. The dense (tenu- 
ous) plasma is resolved by 100 (50) particles per cell per 
species. We employ periodic boundary conditions in all 
directions and no particles are introduced into the sim- 
ulation after t = 0. The plasma clouds detach instantly 
from the boundary at x = ±L x /2 and they will interpen- 
etrate at x = 0. We will discuss in what follows the time 
evolution as shown in the movies followed by the plasma 
state at the simulation's end at the time t = 180. 

In the movies we show the whole time integrated evo- 
lution of the relevant parameters from t = lw~l on to 
t = lS0(jj~] on . Movie 1 shows the separate evolution of 
the field components B x1 B yi B z . Movie 2 shows the evo- 
lution of J z in the upper panel and log | J x + iJ y \ 2 in the 
lower panel. In Movie 1 we see at early times the field 
is planar. The fields have a dipolar structure. As time 
elapses it becomes more and more nonplanar, filamented 
and fragmentary. The B x which is initially uniform as- 
sumes a fragmented morphology, with gradually increas- 
ing structural lengthscales. At early times the forming 
shock accelerates before reaching a steady speed. Striped 
magnetic fields - tell-tale signs of the filamentation insta- 
bility in the B x and B z movies before t < 43u;~* on . Grad- 
ually the magnetic fields become dominated by large- 
scale structures which represent the extrema. B Z1 a 
tracer of the plane current takes on the almost circu- 
lar structure characteristic of the ring current. In Movie 
2 we see the gradual breakup of the initial current sheet 
and the formation of a current ring. The sheet initially 
breaks up into smaller rings at t ~ S2oj~l on which grad- 
ually increase in volume and merge (t ~ 113cj~* on ) to 
form the final large rings, limited only by the box size. 
Incoming current filaments are deflected from the current 
sheet and shear off in the negative y direction, rolling up 
into vortex structures which later merge and distort the 
current sheet. 

Figure [2] displays the phase space distributions / e (#, T) 
of the electrons and fi(x,p x ) of the ions, which have been 
integrated over all y. The f e shows us the energies the 
electrons reach at the shock, while the merging of the 
ions of both clouds in fi(x,p x ) is a necessary condition 
for a plasma shock. A downstream region is visible in 
fi for 55 < x < 70, which separates a reverse shock 
from the forward shock moving to increasing x. The re- 
verse shock is weaker than the forward one, due to the 
asymmetric plasma densities that yield a net flow of the 
downstream region. In what follows we discuss only the 
forward shock. Electrons reach in the cloud overlap layer 
values of V « 100 in significant numbers and the fastest 
ones reach V w 250. An ion with the mass rrii = 250 that 
moves with has in the simulation frame the kinetic 
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FIG. 2: The electron phase space density in the x — Y plane 
(a) and the ion phase space density in the x — p x plane (b) : 
The r > 200 reached by some electrons implies that their 
kinetic energy is comparable to that of the ions. The ion 
distribution reveals a forming downstream region 55 < x < 
70 that separates a reverse shock from the forward shock, 
which moves to increasing x. The ions of the tenuous cloud 
are reflected back into the upstream region within 70 < x < 
80. This interval coincides with that of the strongest electron 
acceleration. The color scale is 10-logarithmic. 



energy of an electron with T ^ 70. A shock transition 
layer with 70 < x < 80 reflects the incoming ions of the 
tenuous cloud and accelerates its electrons, as observed 
at the Earth's perpendicular bow shock |20| . 

Figure [3] depicts B, which has grown in amplitude 
by an order of magnitude compared to |Bo| ~ 16. All 
magnetic components reveal elongated structures in the 
downstream region 55 < x < 70, which have a thickness 
of up to one ion skin depth. This magnetic field is tied 
to the current channels driven by the filamentation in- 
stability and the Bo deflects the particles and thus the 
currents in all directions, causing them to form three- 
dimensional distributions. A massive coherent magnetic 
structure is present in the shock transition layer with 
70 < x < 80. The B x and B y form a vortex, which ro- 
tates counterclockwise in the x — y plane and is centered 
at (x,y) « (74,3), where its modulus \B X + iB y \ « 0. 
The vortex fills up the full ^-interval of the simulation 
box at this time. The B z has its peak at the center of 
the magnetic vortex and its amplitude decreases, radially 
from the centre (x,y) ~ (74,3). 

The magnetic field shows a clear subdivision between 
its components in the simulation plane and the compo- 
nent orthogonal to it, which should hold also for the cur- 
rents. Figure S] shows the currents in and out of the sim- 



3 



(a)B x ,t= 180GW 1 



\ \ 


k 


i 

1 


50 55 60 65 

XA, ions 


70 75 80 


i.O -51.5 -18.9 13.6 
(b)B y ,t= 1 


46.1 78.6 11 

80GW 1 




FIG. 3: The magnetic field components B x (a), B y (b) and 
B z (c): All components evidence structures on ion skin depth 
scales within the 55 < x < 70. The B x and B y reveal a mag- 
netic vortex with a counterclockwise sense of rotation within 



70 < x < 77 and with its center with \B X + iB y 



at 



(x,y) ~ (74,3). A further quasi-planar structure in B y is lo- 
cated within 77 < x < 80. The B z is strongest in the interior 
of the vortex in the x — y plane and its perimeter coincides 
with the magnetic vortex. 



ulation plane. In agreement with the magnetic structure 
in Fig. [3j the currents in Fig. [4] are filamentary in the 
forming downstream region 55 < x < 70 and in the inter- 
val 76 < x < 80, while they look fundamentally different 
in the shock transition layer in 70 < x < 76. Here J x 
and J y form a vortex centered at (x, y) ~ (74,3), while 
J z > is strong in the ellipsoidal interval encircled by the 
vortex. Current channels almost parallel to y are located 
within 76 < x < 80, which give rise to the quasi-planar 
By at this position in Fig. [3jb). These current striations 
are what is left of the current sheet, which developed 
due to the different magnetic deflection of the upstream 
electrons and ions by the perpendicular magnetic field 
component of the shock. 

Figures [3] and 2] reveal why the magnetic field struc- 
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FIG. 4: The currents log 10 (J^ + J;) (a) and J z (b): The 
distribution in (a) reveals an almost circular vortex for 71 < 
x < 76 and —0.5 < y < 5 (periodic wrap-around). This 
vortex encircles an interval with a strong J z > 0. The currents 
show a correlation indicating that it flows in the plane on the 
perimeter but is increasingly deflected along z as we go to the 
center. A current system is present in all components in the 
downstream region 55 < x < 70 and in 76 < x < 80. 



ture in the shock transition layer appears to be stable. 
The current vortex in the x — y plane is responsible for 
B z . This current is almost aligned with the magnetic 
field in the simulation plane in a region with B z w 0, 
as Figs. [3fb) and HJa) show. The current vortex thus 
experiences negligible magnetic force. Its sense of rota- 
tion is counter-clockwise and it gives an interior magnetic 
field that points in the positive z-direction. The B z out- 
side the current vortex is negligible and this system is 
analogous to an infinitely long coil with a symmetry axis 
parallel to z. The current J 2 , which is responsible for 
the vortex formed by B x and B yi also flows along B, 
because \B x J r iB y \ w in the centre of the magnetic vor- 
tex. However, the magnetic field configuration is prob- 
ably only approximately force- free, because the current 
vortex in Fig. H^a) spirals into the interval with B z ^ 
thereby coupling the current components. A force- free 
magnetic field must fullfill V x B = aB with a constant 
a [2l|, [22| , which is achieved here by the combination of a 
solenoidal component and of a ring component of B with 
a comparable strength. 

Discussion: Our simulation demonstrates the growth 
of the filamentation instability, in spite of the initial con- 
ditions that reduce its growth rate relative to the com- 
peting instabilities [l6j|. The filamentation instability is, 
however, not the main source of the magnetic field within 
the shock transition layer. A current sheet develops at 
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the front of the dense cloud, because the electrons and 
ions are not equally deflected by the perpendicular com- 
ponent of B. This current sheet boosts the component 
of B orthogonal to the flow direction to an amplitude 
& lOO^o^, resulting in a magnetic energy density that 
is comparable to the initial ion kinetic energy density. 
This current sheet is unstable and a magnetic structure 
develops, which is composed of solenoidal and ring com- 
ponents. Such a topology can, in principle, be force-free. 
However, the current loop in the simulation plane is not 
closed and it is also elliptical. In fact, we cannot ex- 
pect this magnetic field structure to be force-free. The 
ram pressure of the inflowing upstream plasma, the elec- 
tron acceleration to extreme speeds and the ion reflection 
in the shock transition layer exert a force on this struc- 
ture and they are probably responsible for its compres- 
sion along x into an ellipse. The spiral structure of the 
current in the simulation plane may be a remnant of the 
growth of the structure and thus a transient effect. 

Clearly our simulation results are affected by the pe- 
riodic boundary conditions in the y direction and by the 
spatial 2D geometry. The 2D geometry assumes that 
the magnetic structure is an infinitely long cylinder with 
the symmetry axis along z. A force- free 3D topology 
is obtained if the cylinder is bent to form a closed loop 
known as a spheromak, which is composed of poloidal 
and toroidal magnetic fields [U, [22| . The structure we 
observe in our simulation would then be a cross-section 



of the torus. The obvious effect of the periodic boundary 
conditions along y is to halt the growth of the magnetic 
structure, by which a steady state is reached (see the 
movies). This steady state is actually quite remarkable, 
considering the powerful particle acceleration processes 
taking place at this location. A wider y interval will pre- 
sumably result in a further growth of the magnetic struc- 
ture. Spheromaks can be quite large and they have been 
invoked to explain large scale structures in the solar wind 
triggered by coronal mass ejections [23|. The poloidal 
magnetic field of the spheromak, which corresponds to 
B z in our simulation, is strong and homogeneous within 
the torus' perimeter. 

Our simulations suggest that spheromaks, or similar 
magnetic loops, that form within the current sheet of 
an internal GRB shock may thus provide the means to 
grow stable structures with a strong coherent magnetic 
field from kinetic scales to MHD scales. The current 
sheet, which is the source of this structure, does require 
an ambient magnetic field with a component orthogonal 
to the flow velocity vector and ions. 
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